VERIFICATION OF A NEW NON-LINEAR IV-EXPONENT: 
SIMULATION OF THE 2D COULOMB GAS WITH LANGEVIN DYNAMICS. 
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0\ . ABSTRACT 

' It has recently been suggestedi from scaling arguments that the non-linear IV- 

. exponent a, V (X I'^, for a two-dimensional superconductor is different from the 

1^ ' exponent originally suggested by Ambegaokar et al.(AHNS)u. The relation between 

pLj . the new and the old exponent is a = IttAHNS — 3. The new scaling behaviour 

I is linked to the logarithmic vortex interaction and the long range time tail which 

CSJ . this gives rise to. Consequently one may expect that the scaling behavior is generic 

for models which have these basic features. The simplest model of this type is the 
two-dimensional Coulomb gas model with Langevin dynamics. We here explicitly 
verify, through computer simulations, that the /^-characteristics of this model 
. indeed scales according to the new scaling exponent a. 

o : 

^ ■ The transition from the resistive state to the superconducting state is for two 

ly-^ , dimensional (2D) superconductors of the Kosterhtz-Thouless (KT) type.! In the su- 

0> \ perconducting state the current- voltage {IV) characteristics is non-linear, V ^ 1°^, 

' where the exponent a is larger than one. An expression for a was obtained by Am- 

O . begaokar et al (AHNS) in ref.lil. However, it has recently been suggested that the 

I ! exponent a is linked to a critical dynamics of the vortex fluctuations in the low tem- 

p I perature phase and that, as a consequence, it is different from the AHNS-value.H The 

^ ■ relation between the new and the old value of the exponent is a = 2aAHNS — 3.0 

•• • The correctness of the new exponent was corroborated by simulations of the 2D RSJ- 

• w-{ ', model.0 Since the new exponent is linked to a critical dynamics and a critical scaling 

^ I argument, it is interesting to know what ingredients are necessary to get this critical 

■ dynamics. In the present paper we find a candidate for the simplest generic model 

that has the required critical property^. 

It is well-known that the KT-transition is due to unbinding of thermally created 
vortex-antivortex pairs and that the voltage is due to flux flow of free vortices. H The 
interaction between the vortices is logarithmic with distance. Our hypothesis is that 
the critical dynamics is linked to the logarithmic vortex interaction. The Coulomb 
gas interaction is logarithmic in two dimensions and consequently the thermally cre- 
ated vortices and antivortices can be described as a two dimensional gas of Coulomb 
charges.!! The vortices are quantized and can have vorticity s = ±1, which corre- 
sponds to positive and negative Coulomb gas charges. The total vorticity is always 
zero which means that the Coulomb gas always has equally many positive and neg- 
ative charges. The vortices of a 2D superconductor obey a particular dynamics.El 
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Consequently the vortex part of a 2D superconductor can be modelled by a Coulomb 
gas with a particular dynamics for the Coulomb gas charges. However, if the critical 
dynamics is strongly linked to the logarithmic interaction, one may suspect that it is 
relative insensitive to the dynamical equation of the moving particles. This suggests 
that a 2D Coulomb gas with a simple dynamics would fall in the class of models 
which have the required critical property. In the present paper we show, through 
computer simulations, that a 2D Coulomb gas with conventional Langevin dynamics 
indeed belongs to this class of models. 

The equation of motion for the 2D Coulomb gas with Langevin dynamics is given 

by 

^ = fF,,W + „W (1) 

where r is the position and Ftot is the total force acting on it due to all the other 
particles as well as any externally imposed force, D is the diffusion constant, T is 
the temperature (unit system such that the Boltzmann constant fc^ = 1), and r/ is a 
random force obeying 

{r]"{t)r]^{t')) = 2D6^p6{t-t') (2) 

where a and /3 denote the Cartesian components. The charges are circular disks with 
diameter ro such that the force acting between two particles i and j with charges 
Si and Sj respectively (units such that the charge is s = ±1) and separated by the 
distance r is given by 

Fij = SiSj{---Ki{r/ro)) (3) 
r ro 

where Ki is a modified Bessel function of order 1. Note that the force between two 
particles vanishes for r = 0. This means that the charge distribution of a particle is 
soft, which is in accordance with the precise vortex-Coulomb gas particle analogy.^ 
In our simulations we choose a 2D box of length L with periodic boundary conditions 
and a constant particle density n. The results presented here are basically for n = 
5 X IO^^Tq^ and L/tq = 226 ( i.e. the total number of particles is = 256). The 
simulations are done by discretizing time with a time spacing At and introducing a 
Gaussian random noise T]{t) acting independently on each particle at each time step. 
The equation of motion is then solved on the computer by using a standard Euler 
integration method.cl The largest time sequences used in the present simulations are 
given by t/ At ^ 5 X 10^. The At chosen has to be small enough to ensure that the 
equation of motion is correctly solved yet as large as possible in order to achieve as 
large time sequencies as possible.a 

We are calculating the average charge current Ip as a function of an external force 
F^xt- The average charge current is given by 

^P = (E^^^) = ^(EFl:kt)) (4) 

i=l "-^ ^ i=l 

where the sum is over all particles and F^o^ is the total force acting on the particle 
i and the brackets () denote a time average. The voltage \^ in a 2D superconductor 



is due to the flux flow which means that V is proportional to Ip. The force acting 
on a vortex is given by the Lorentz force and is consequently proportional to the 
current. This means that the externally imposed current / is proportional to the 
external force Fex- Thus the JV^-characteristics for a 2D superconductor corresponds 
to the Fea;/p-characteristics of the Coulomb gas. Figure 1 shows the obtained F^^Ip- 
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Fig. 1. log-log plot of the Fexlp characteristics (T = 0.12, 0.14, 0.16, 0.18, 0.23, 0.26, from bottom 
to top, respectively). The pluses represent the data and the dashed lines are fits to the form 
AKniBFex) + C as described in text. 



characteristics for a sequence of temperatures plotted as In Jp against InFe^^. As seen 
in the figure the data fall on straight lines in the limit of small F^x which implies 
that Ip oc F^^ where the slope of the line gives the corresponding exponent a. In 
this way the non-linear F^xlp- exponent a is obtained. As also seen in figure 1 the 
Fexlp-data. can be very well fitted to the functional form AKo{BFex) + ClnFex over 
an appreciable range of F^x, where A = 1 — a, and B and C are constants. This 
implies that Ip = CFexe^p{AKQ{BFex)) and this observation somewhat simplifies the 
determination of the exponent a.0 The obtained exponents a, corresponding to the 
data in figure 1, are shown in figure 2 together with estimated error bars. According 
to the critical scaling theory these exponents a should be related to the static charge 
density correlations.H The precise relation is 

a=^-l (5) 

where e is related to the Fourier transform of the dielectric function e(k). The dielec- 
tric function is in turn related to the charge density correlations by 

i^ = l-^(Mk)An(-k)) (6) 



7.0 



(U 

a 

o 
a. 

X 
(U 



_ 




1 


1 1 

A MP scaling 




\ 




V AHNS 


- 




\ 


• Langevin simulations 


- 






































-^-^-^-^ 








^^^^^ J 

1 



> 3.0 



1.0 

0.10 0.15 0.20 0.25 

Coulomb Gas Temperature, T 

Fig. 2. The non-linear IV-exponent a as a function of temperature (dashed lines are guides to the 
eye). 

where 

N 

An(r)=^s,5(r-ri) (7) 

is the charge density. The Fourier transform of the dielectric function has the leading 
small k dependence^ 

-1- = ^ (H) 
e(k) eP + A-2 ^ ^ 

where A is the screening length. This relation defines e in terms of e(k). The screening 
length A is infinite for temperatures below the KT transition temperature. This means 
that e = e{k = 0) in the low temperature phase. 

In the simulations we calculate An(k, t) for a long time sequence for a given 
fixed value of k. From such a sequence we readily obtain the time average of the 
charge density correlations (Ari(k)An(— k)) and from this l/e(k). The value of e 
is then determined by using equation (^. This determination of e is shown in 
figure 3 and figure 4 gives the values obtained . The scaling exponent a = 1/Te — 1 
is now calculated from these values of e and is compared to the actual exponent a 
determined directly from our /pFg^-simulations in figure 2. As seen in figure 2, the 
agreement with the scaling theory prediction is excellent. According to the AHNS- 
theoryS the exponent a should instead be given by aAHNS = l/2Te -|-l = a/2-|-3/2. 
This prediction is also shown in figure 2. As is apparent from figure 2, the AHNS- 
prediction does not agree with the simulated exponents. Thus figure 2 clearly verifies 
that the new non-linear IV exponentil gives a correct description of the 2D Coulomb 
gas with Langevin dynamics and that the AHNS-predictionS is not correct for this 
model. The same conclusion was obtained in case of the 2D RS J-model.i The present 
simulations suggests that the new scaling exponent is quite robust and that particles 
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Fig. 3. The dielectric function as a function of |k| for T=0. 12,0. 14,0. 16,0. 18,0. 23,0.26, from top to 
bottom, respectively, and the corresponding fits used to obtain e. 
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Fig. 4. 1/e obtained by fitting equation (||) to the data in figure | 



with a logarithmic interaction seems to be the most essential feature for obtaining 
this new exponent. 

According to the scaling theoryi, the new critical exponent is linked to a long 
time tail in the charge density correlations. We have explicitly caculated the long 
time behaviour of the charge density correlations function 

^(k,t) = (n(k,t)n(-k,0)) (9) 

From a practical point of view it is easier to converge the average for large values of 
t if one uses a /c- value somewhat larger than the smallest possible k = 2tt/L. Figure 
5 shows the result for a temperature (T = 0.18) somewhat below the KT transition 
{Tkt ~ 0.22). In this case a resonable convergence was achieved for the k- vector 
k = ^(4,4). The result is plotted as tg{\i,t) versus t. As seen in the figure tg{^^,t) 
is constant to very good approximation for large t which means that ^(k, t) oc 
The result given in figure 5 shows that the time correlations for a small finite k fall 
to good approximation like 1/t over a very large t-interval. This result suggests that 



g{li = 0, t) has a true 1/t tail for large t. 

In summary the results displayed in figures 1-4 demonstrate that the 2D Coulomb gas 
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Fig. 5. The charge density correlation function as a function of t. The data clearly shows that 
the charge density correlations decay as 1/t over a large time interval. Here t is in units of 
r'^/D {At = 0.08 rl/D). 

with Langevin dynamics has a non-linear /\^-exponent a which is distinctly different 
from the AHNS-predictionl^ and is instead given by the scaling theory by Minnhagen 
at al.a. In addition the result given in figure 5 is indeed consistent with a 1/t-tail in 
the charge density correlations, as required by the scaling theory. We expect these 
findings to be quite general and so should apply to a 2D superconductor and a 2D 
Josephson coupled array. 
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